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Abstract 

We compute the electroweak corrections due to the light fermions to the pro- 
duction cross section a{gg — * H) and to the partial decay widths V (H ^77) and 
r [H — > g g). We present analytic results for these corrections that are expressed in 
terms of Generalized Harmonic Polylogarithms. We find that for the gluon fusion 
production cross section and for the decay width V (H — > gg) the corrections are 
large in the Higgs mass region below 160 GeV where they reach up to 9% of the 
lowest order term. For the decay width V (H — > 77) the corrections for Higgs mass 
above 160 GeV can reach —10% of the lowest order term. 
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One of the great successes of the experimental program carried out at LEP has been 
to put a firm lower bound on the Higgs mass, m H > 114 GeV £Q, and at the same time, 
together with the information coming from SLD, to give a strong indirect evidence that the 
Higgs boson, the still missing particle of the Standard Model (SM), should be relatively 
light with a high probability for its mass to be below 200 GeV. The search for the Higgs 
boson is one of the main objective of the Tevatron and the future Large Hadron Collider 
(LHC), that are supposed to span all the Higgs mass regions up to 1 TeV. At hadron 
colliders the main Higgs production mechanism is the gluon fusion [2|, a process whose 
knowledge is fundamental in order to put limits on the Higgs mass or, in case the Higgs 
is discovered, to compare the measured cross section with the SM result. Concerning 
the Higgs decay channels, it is quite difficult for an hadron collider to access part of the 
mass range favored by the LEP results, the so called intermediate Higgs mass region 
114 < m H < 160 GeV, because of the large QCD background to the dominant modes. In 
this region the rare decay H — > 77 is the most interesting alternative to the usual decay 
channels. 

It is then natural to try to have a theoretical prediction of the gluon fusion production 
cross section as well as this decay mode, as precise as possible. In the recent years most of 
the attention to these processes has been devoted to the calculation of the QCD corrections, 
but we have now reached a stage such that also the electroweak effects can be interesting. 
In this report we make a first step in the calculation of the electroweak corrections by 
presenting the analytic result for the two-loop contributions induced by the light (assumed 
massless) fermions. This subset of corrections, finite and gauge invariant, is actually 
interesting for two reasons: i) it could be numerically not irrelevant because one has to 
sum over the generations; ii) because of the presence in the diagrams of two mass scales, 
m H and m w , that in the most physically interesting Higgs mass range are relatively close, 
and of cuts connected to the zero mass fermions, this subset cannot be computed with the 
standard heavy mass or momentum expansion method but requires a complete calculation, 
making it a good probe of new techniques for the computations of two-loop integrals. 

We start analyzing the production cross section a(gg — > H). Because of its importance 
this process has been investigated with great accuracy in the recent years and it is now 
known at the next-to- next-to-leading order in QCD jSJIUE]. A recent discussion [S] on 
the residual theoretical uncertainty from perturbative QCD contributions estimate it to 
be below 10% for m H < 200 GeV. Electroweak corrections to this production mechanism 
were only considered in the large-mj limit and found to give a very small effect, below 1% 



The Higgs boson, carrying no color, has no tree-level coupling to gluons; therefore this 
process proceeds via loops. At the partonic level the cross section, not corrected by QCD 
effects, can be written as: 



where the lowest order one-loop contribution is only due to the top quark and is given by: 
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Figure 1: Two-loop topologies involving light fermions contributing to gg — > H. 
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At one-loop the contribution of light fermions is suppressed by their coupling to the Higgs 
and it is completely negligible. 

At the two-loop level the light fermions can contribute to this process because their 
Higgs coupling suppression can be avoided by coupling them to the W or Z bosons that 
can directly couple to the Higgs particle. Indeed the light fermions contribute to the Higgs 
production through the topologies (a) and (b) depicted in Fig. HI with the W or Z boson 
exchanged in the loops together with quarks. 

The general structure of the amplitude for the production of a Higgs particle, in the 
fusion process of two gluons of polarization vectors 6^(51) and e u (q 2 ), can be written as: 

= q? q\ T X + q» q v 2 T 2 + q» 2 T 3 + q\ q» T 4 + (q, ■ q 2 ) gT T 5 + <?>"" q lp q 2(7 T 6 , (4) 

with gauge invariance requiring that T\ = T 2 = and T4 = — T5 while T3 does not con- 
tribute for on-shell gluons and Tq vanishes exactly. The relevant form factors T4 and T5 
are extracted from the one-particle irreducible (1PI) diagrams of Fig. ^with the use of 
standard projectors. The evaluation of the 1PI part of the form factors is performed in 
few steps. First, the scalar amplitudes are reduced to a set of linearly independent ones. 
This step is achieved in two different ways jS]: i) a reduction in terms of scalar integrals 
containing only denominators is obtained with the introduction of a fictitious auxiliary 
planar double box diagram 9\, with final momenta equal to the initial ones, that contains 
an "extra" propagator with respect to the diagrams of Fig. HJ such that the seven possi- 
ble invariants that can appear in the numerator of a scalar amplitude can be simplified 
against the seven propagators of the auxiliary diagram, ii) By shifting the loop momenta 
in an appropriate way we write the scalar amplitudes in terms of integrals with a number 
of propagators equal or smaller than the original one, having a set of independent scalar 
products in the numerators [10j . Next, the independent scalar amplitudes obtained in the 
two reduction methods are expressed through the use of integration-by-parts identities 
[TT] as a linear combination of the same finite set of basic integrals, called master integrals 
(Mis). The subsequent evaluation of the Mis is performed employing the technique of the 

1 All the analytic continuations are obtained with the replacement x — * x — i e 
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differential equations in the external kinematical variables ^2] • Finally, the coefficients of 
the Mis Laurent expansion in e = (4— n)/2, where n is the dimension of the space-time, are 
expressed in closed analytic form. The 1PI diagrams of Fig. ^ show thresholds in the ex- 
ternal Higgs momentum, q 2 , at q 2 = 0, m^m^. Because of the presence of the q 2 = 
and q 2 = 4 w? w thresholds our analytic result cannot be expressed only in terms of Har- 
monic Polylogarithms (HPLs) JHIE], a generalization of Nielsen's polylogarythms, but 
requires a recently introduced extension of the HPLs, the so-called Generalized Harmonic 
Polylogarithms (GHPLs) [T3]. We explicitly checked that our analytic result satisfies all 
gauge-invariance requests. 

The total light fermion contribution to the partonic production cross section <r(gg — > 
H) is obtained summing over the first two generations of quarks both the diagrams with 
the W and those with Z bosons. For the bottom quark instead, we can include only the 
contribution due to the Z boson, which is by itself finite and gauge invariant. We find, in 
the same normalization of Eq. (0) and in units a/(2irs 2 )(m 2 v /m 2 v 
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where w H = mr w jrr? H , z H = m 2 /m H , s = sin 0w, c = 1 — s and 
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Eq. © is expressed in terms of ordinary HPLs, the H functions with only indices 0, — 1, 
and of GHPLs, H functions where also the indices — r and —4 are present. Their explicit 
expressions are presented in Appendix A. 

In Fig. El we plot the relative corrections 5, induced by the two-loop light fermion 
contribution, to the Higgs production cross section o = oo(l + S), where <To is the lowest 
order result. The correction raises till the opening of the 2W threshold in the topologies 
(a, b) of Fig. [H reaching up to 9% of the one-loop result. As soon as the threshold is 
passed the correction has a sharp decrease with only a small bump at the opening of the 
2 Z threshold, and for values of m H larger than 200 GeV remains always smaller than 
2%. The large (for an electroweak correction) result found in the intermediate Higgs mass 
region can be explained noticing that the topologies (a, b) of Fig. [I] with their thresholds 
at 2 W and 2 Z are of a new type with respect to the one-loop result that is described by 
a top loop triangle diagram. We notice that, in the region below the 2W threshold, the 
size of the electroweak light fermion corrections is comparable to the accuracy reached in 
the knowledge of the QCD corrections to this production cross section jH]. 

The techniques used to compute the two-loop light fermion corrections to the gluon 
fusion amplitude can be applied to study also the decay processes H — > 77, gg. The 
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Figure 2: Relative corrections to the production cross section a(gg — ► if) and to the 
decay width r (H — > gfgr). 



result for the latter is very easy to derive because one can relate the partial decay width 
r(# -> flfflf) to a (gg ^ H) via: 

Y{H^gg) = ^a(gg^H). (7) 

Thus our result for the light fermion corrections to the production cross section can be 
translated into a result for the corrections to the partial decay width T(H — > gg) and Fig- EI 
can also be read as the relative corrections to this decay mode. It should also be noticed 
that the QCD corrections to this decay are extremely large jTHJH]. As an example, in the 
intermediate Higgs mass region they shift the decay width by about 60-70% upward. 

We consider now the decay H — > 7 7 that is of particular interest for the Higgs searches 
at hadron colliders thanks to its clean signature. The Higgs boson, being a neutral particle, 
has no tree-level coupling to photons, therefore also this decay can proceed only via loops. 
The corresponding partial decay width can be written as: 

r #^77 = _ \t 2 , 8 

where the lowest order one-loop contribution is due to W boson and top loops, and is 
given by [T7] : 

F U = rZ + ^N e F?, (9) 
J* = 2(l + 6w H )-12w H (l-2w H )H[-r,-r;-—) , (10) 



t} 1 = gu, (11) 
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Figure 3: Two-loop topologies involving light fermions contributing to H —>■ 77. 

with N c the color factor. 

The W and top one-loop contributions interfere destructively. They approach constant 
values {T^ 7, T} 1 — > —4/3) for mass of the particle inside the loop much heavier than 
m H , with the W loops always providing the dominant part for Higgs mass up to m H ~ 600 
GeV, where the top contribution becomes comparable. 

Two-loop investigations of this decay channel have been mainly devoted to QCD cor- 
rections to the one-loop top contribution |18j . Concerning two- loop electroweak corrections 
only the 0{G tl m1) term [TH] and the 0{G^m 2 H ) one [20] were studied. 

Differently from the one- loop level at two- loop the light fermions can contribute to 
this partial decay width through all diagrams whose topologies are depicted in Fig. El 
The choice of an appropriate gauge-fixing procedure can somewhat reduce the number of 
diagrams to be computed. To achieve this we find it convenient to employ the Background 
Field Method (BFM) quantization framework. The BFM is a technique for quantizing 
gauge theories j2H 122] that avoids the complete explicit breaking of the gauge symmetry. 
One of the salient features of this approach is that all fields are splitted in two components: 
a classical background field V and a quantum field V that appears only in the loops. The 
gauge-fixing procedure is achieved through a non linear term in the fields that breaks 
the gauge invariance only of the quantum part of the lagrangian, preserving the gauge 
symmetry of the effective action with respect to the background fields. Thus, the BFM 
gauge-fixing function induces a modification of some of the vertices in which background 
fields are present, such that, in the BFM Feynman gauge (£q = 1, £q being the gauge 
parameter), the vertex coupling a background (external) photon to a quantum W boson 
and its unphysical scalar counterpart is absent 2 . With this choice of gauge only the physical 
W and Z bosons are exchanged in Fig. 01 In particular, in the topologies (a, b, c) the 
internal wavy line can only be a W boson while in the topologies (d, e) it can represent 
both a W and a Z boson. The evaluation of the diagrams of Fig. |3] follows the line already 

2 A general discussion on the application of the BFM to the SM with the complete set of BFM Feynman 
rules can be found in Ref. while a specific two- loop calculation in this framework is presented in Ref . [21] ■ 
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discussed for the gluon fusion cross section. Concerning the counterterm diagrams, which 
are not depicted in Fig. their evaluation does not present any particular difficulty. We 
just notice that the contributions coming from the 0(e) part of the one-loop diagrams 
cancel out in the total result. The renormalization is actually straightforward. Because 
we are in the presence of two real photons their coupling is given by a = 1/137.036 while 
the coupling of the Higgs to the vector boson is more appropriately described in terms of 
Gy,. Correspondingly, we employ the on-shell definition for the W boson mass. 

The full light fermion contribution to the Higgs decay into two photons is obtained 
summing all the diagrams over the 3 lepton families plus the two first generations of quarks. 
Concerning the third generation, the bottom contribution to topologies (d, e) mediated 
by the Z boson, that is finite and gauge- invariant, can also be included. The light fermion 
contribution to the partial width H —>■ 77, in the same normalization of Eq. (jHJ), reads in 
units a/(2Tis 2 )(m 2 v /m 2 H ): 

3f f = 2tf c A 2 [-2/9,^] + 3A 2 [0,i<; ff ] + 

+4 (l-2s 2 + As 4 ]A 1 [z H ] , (12) 



where 

A 2 [q,x) = -8(l + q)+A(l+q)(l-x)H (-1;--) -2(l + 2qx)H [0,-l;- 1 



— (5 — 12 x) H ( — r, — r; j — 6 (1 + q — 3x — 2 qx) H \—r, — r, — 1; — 

3 \ xj \ x 

2x) Hi 0,-r, -r:-- ] +(l-3x)H\ 0,0,-1:-- 



+2 (1 + 2?) 

-Vl -4x |2 (1 + 2 q) H f-r; ~ j - 6 q (1 - 2x) H (-4, -r, -1; 



+4q (1 - 2x) 



H I -r, 0, -1; 1 + H ( -r, -r, -r; — 



6(1 - 2x) 2 „ ( I 



H —r, -1; — , (13) 



VI -4a; \ x 

and Ai has been given in Eq. (|5jl. 

Eq. (fT3*|l shows an unphysical singularity at x — 7, i.e. at the opening of the 2W 
threshold, connected to the appearance in its last line of a term with a square root in 
the denominator that exactly at the threshold becomes formally infinite. This unphysical 
infinity is actually a signal that our first order treatment of the W propagator in topologies 
(a, b) of Fig. El is inadequate in the 2 W threshold region. To obtain a finite result also in 
this region we regulate the singularity by performing the replacement m w — > m w — iT w /2 
in the square root and we check the dependence of the final result on the regulator. 

In Fig. El we plot separately the relative corrections, 5, to the decay width T = Tq(1 + 
5) (To oc l^ 1 '] 2 ) induced by the two-loop light fermion contribution and by the QCD 
corrections to the one-loop top term as well as the sum of the two contributions. The 
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Figure 4: Relative corrections to the decay width H —>■ 77. 



figure is obtained by regulating the singularity at x = 7 with Y w /m w = 2.5 ■ 10~ 2 . We 
check the dependence of the result on the regulator by varying Y w /m w down to 2.5 ■ 1CT 9 . 
Above the threshold already at m H = 165 GeV the result changes by less than 5% and 
becomes insensitive to the size of the regulator at m H = 170 GeV. Below the threshold, 
down to m H = 150 GeV the variation is larger. However, for a Higgs mass value around 
150 GeV the light fermion correction is small and therefore it is not surprising that the 
relative variations are large. It is clear that in the region 150 < m H < 170 GeV our result 
should be taken with some caution. 

As it can be seen from Fig. 0] the QCD and light fermion contributions for m H < 140 
GeV are of comparable size but with opposite sign so that the total correction is below 1%. 
As an example, for m H = 120 GeV the light fermion contribution amounts to ~(— 1.5%) 
of the total width while the QCD one gives ~(+1.8%). For values of m H above 140 GeV 
the light fermion correction becomes small until the region close to the opening of the 
2 W threshold is reached. Crossing the threshold the correction becomes quite large, in 
absolute value, but as soon as the 2 W threshold region is passed it starts to decrease 
again. However, for m H around 200 GeV, we find still a ~5% correction while the QCD 
contribution is relatively small. The visible bump in the figure around 180 GeV is actually 
connected with the opening of the 2 Z threshold in the topologies (d, e) of Fig. EJ 

In conclusion, we have computed the two-loop electroweak corrections induced by the 
light fermions to the Higgs gluon fusion production cross section o (gg — > H) and to the 
partial decay widths r (H — > 7 7) and r (H — > g g). For the gluon fusion production cross 
section we have found that the light fermion corrections are quite large in the region 
below the 2 W threshold reaching up to 9% of the lowest order result. In the case of 
the Higgs decay to two photons, the corrections are large above the 2 W threshold, while 
below it they are quite small a part the region m H < 140 GeV where they are larger 

7 




2-loop full — 
2-loopQCD - 
2-loop EW 



than 1% cancelling part of the QCD corrections giving a total correction below 1%. It 
should be recalled that, both for the production and the decay, the light fermions are only 
a small part of the electroweak corrections. A complete calculation of the electroweak 
corrections to cr(gg — > H) seems to be relevant to assure the accuracy of the knowledge of 
this important production cross section. 

One us (G.D.) wants to thank F. Maltoni for interesting discussions. This work was 
partially supported by the European Community's Human Potential Programme under 
contract HPRN-CT-2000-00149 (Physics at Colliders). 



A Generalized Harmonic Poly logarithms 

In this Appendix, we outline the definitions and properties of the HPLs and GHPLs and 
give the explicit expressions of those entering in Eqs. (|6I13|) . 

The idea lying behind the introduction of HPLs is to express a given integral coming 
from the calculation of a Feynman diagram in a unique and non-redundant way as a linear 
combination of a minimal set of independent transcendental functions. These functions 
are expressed as repeated integrations over a starting set of basis functions and this set 
depends strongly on the problem one has to solve, being connected directly to the threshold 
structure of the diagrams under consideration. In the case in which only thresholds at 
and 4m 2 or at and m 2 occur, the original set of HPLs introduced in is sufficient to 
cover all the possible integrals. However, for different structure of the thresholds of the 
Feynman diagrams involved in the calculation, as it is the case at hand, one is forced to 
enlarge the set of starting basis functions. 

In our case, we enlarge the original basis functions with indices {—1, 0, 1}: 

111 

g(-l-x) = — — , g(0;x) = -, g(l;x) = - ; (14) 

1 + X x 1 — X 

that give a "weight- 1" HPLs: 

H(-l;x) = -log(l -x) , H (0; x) = log x , H(l; x) = log(l + re) , (15) 

with the following two basis functions [To] : 

1 1 

ff(-4; x) = —— , g(-r; x) = , (16) 

4 + x y/ x {l + ar) 

with the corresponding weight- 1 GHPLs: 

H(-4;x) = log(4 + x) - 2 log 2, H(-r; x) = - logz , (17) 

where: 

z= Vx+Z-Vx 

Vx + 4+y/x 

Integrals with weight > 1 are defined using the recursive definitions: 

i r 

H(0 W ; x) = - log™ x , H(a,w;x)= / dt g (a; t) H (w; t) , (19) 
w Jo 
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where O w = (0, 0, • • ■ , 0) is a vector containing w zeroes, w is a vector with w components 
running over the set {—4, — r, —1, 0, 1} and a takes the values —4, — r, —1, 0, 1. 

It has been shown [T3J [HJ ES] that the set of (G)HPLs constructed in this way satisfies 
useful relation among themselves. As an example, the product of two GHPLs of weight 
Wi and u>2 can be expressed as a linear combination of GHPLs of weight W\ + w 2 . 

Eqs. (|6I13|) contain 11 GHPLs. H(—r, —r;x) is defined in Eq. (JHJ) while H(—l;x) and 
H(—r;x) have already been presented in Eqs. ()15I17|) : the other 8 are: 

H(0,-l;x) = -Li 2 (-x), (20) 
#(0,0,-1; a) = -Li 3 (-x) , (21) 

H(-r,-r,-r;x) = -Uog 3 z , (22) 

6 

H(0,-r,-r;x) = -~ log 3 z + 2 S lj2 (l - z) (23) 

6 

H(-r,-l;x) = / = log(l + t), (24) 

x t(4 + 1) 

H(-rA).-l:.r) = / -=^= C ^ log (1 + q) , (25) 



V"t(4 + t) Jo Q 

H(-r,-r,-l;x) = I JL_ [ . dq log(l + g), (26) 



o v^(4Tt)7o ^/g(4 + g) 
ff(-4,-r,-l;i) = £^/^ log(1 + 9 ). (27 ) 

where Li 2 (x) = — fjf dt log ^~^ ; and Li 3 (x) = J* dt hl2 ^ are the Nielsen's polylogarithms 
and Si j2 (l — z) = \ z dt log ^— ^ is a Spence function. 

These 11 GHPLs are real and analytic in all the half plane x > and they present 
cuts in the plain x < 0. In our case the variable x is defined as x = —s/m 2 , where s 
is the squared cm. energy and m is the mass of the particle running in the loop (W 
or Z bosons). In the physical region s = m 2 H > 0, the variable x is negative, and the 
GHPLs have to be continued analytically. The analytical continuation is done with the 
usual ze-prescription, i.e. giving a small positive imaginary part to s: s + ie. In so doing, 
we have: 

x — > —y — ie , (28) 
where y = s/m 2 > 0, and, for example: 

y/x -> y/-y -ie = -iy/y, (29) 
logx — > log {—y — ie) = logy — in . (30) 

The numerical evaluation of the 11 GHPLs comes straightforwardly from their definition 
as repeated integrals. 
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